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ABSTRACT 

We present deep optical observations of the gravitational lens system CLASS B02 18+357 
from which we derive an estimate for the Hubble Constant (Hq). Extensive radio observations 
using the VLA, MERLIN, the VLBA and VLBI have reduced the degeneracies between H 
and the mass model parameters in this lens to one involving only the position of the radio- 
quiet lensing galaxy with respect to the lensed images. B02 18+357 has an image separation 
of only 334 mas, so optical observations have, up until now, been unable to resolve the lens 
galaxy from the bright lensed images. Using the new Advanced Camera for Surveys, installed 
on the Hubble Space Telescope in 2002, we have obtained deep optical images of the lens 
system and surrounding field. These observations have allowed us to determine the separation 
between the lens galaxy centre and the brightest image, and so estimate Hq. 

We find an optical galaxy position - and hence an Ho value - that varies depending 
on our approach to the spiral arms in B02 18+357. If the most prominent spiral arms are left 
unmasked, we find Hq = 70±5 km s _1 Mpc~ 1 (95% confidence). If the spiral arms are masked 
out we find H = 61 ±7 kms" 1 Mpc _1 (95% confidence). 
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1 INTRODUCTION 

Objects at cosmological redshifts may be multi- 
ply imaged by the action of the gravitational field 
of foreground galaxies. The first such example 
of gravitational lensing was the system 0957+561 
(Walsh, Carswell & Weymann 1979) in which the 
core of a background quasar is split into two im- 
ages 6" apart. Since then approximately 70 cases 
of gravitational lensing by galaxies have been 
found. 1 

Refsdal (1964) pointed out that such multiple- 
image gravitational lens systems could be used to 
measure the Hubble constant, if the background 
source was variable, by measuring time delays be- 
tween variations of the lensed image and inferring 
the difference in path lengths between the corre- 
sponding ray paths. The combination of typical 

A full compilation of known galaxy-mass lens systems is given on the 
CASTLeS website at http://cfa-www.harvard.edu/glensdata 



deflection angles of ~1" around galaxy-mass lens 
systems with typical cosmological distances im- 
plies time delays of the order of weeks, which 
are in principle readily measurable. Time delays 
have been measured for eleven gravitational lenses 
to date: CLASS B0218+357 (Biggs et al. 1999; 
Cohen et al. 2000), RXJ 0911.4+0551 (Hjorth 
et al. 2002), 0957+561 (Kundic et al. 1997; Os- 
coz et al. 2001), PG 1115+080 (Schechter et al. 
1997) CLASS B 1422+231 (Patnaik & Narasimha 
2001), SBS 1520+530 (Burud et al. 2002a), 
CLASS B 1600+434 (Koopmans et al. 2000; Bu- 
rud et al. 2002b), CLASS B 1608+656 (Fassnacht 
et al. 1999; Fassnacht et al. 2002), PKS 1830-21 1 
(Lovell et al. 1998), HE 2149-2745 (Burud et al. 
2002b) and HE 1 104-1805 (Ofek & Maoz 2003). 
In principle, given a suitable variable source, the 
accuracy of the time delay obtained can be bet- 
ter than 5%. This has already been achieved in 
some cases (eg. 0957+561) and there is no doubt 
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that diligent future campaigns will further improve 
accuracy and also produce time delays for more 
gravitational lens systems. 

Gravitational lenses provide an excellent 
prospect of a one-step determination of H on cos- 
mological scales. The major problem is that, in 
addition to the time delay, a mass model for the 
lensing galaxy is required in order to determine 
the shape of the gravitational potential. The model 
is needed to convert the time delays into angu- 
lar diameter distances for the lens and source. In 
double-image lens systems in which the individ- 
ual images are unresolved this is a serious problem 
as the number of constraints on the mass model 
(lensed image positions and fluxes) allows no de- 
grees of freedom after the most basic parameters 
characteristic of the system (source position and 
flux together with galaxy mass, ellipticity and po- 
sition angle) have been fitted. In four-image sys- 
tems the extra constraints provide assistance, and 
in a few cases, such as the ten-image lens system 
CLASS B 1933+503 (Sykes et al. 1998) more de- 
tailed constraints on the galaxy mass model are ex- 
ploited (Cohn etal. 2001). 

There are two further systematic and poten- 
tially very serious problems. The first is that the 
radial mass profile of the lens is almost com- 
pletely degenerate with the time delay, and hence 
H (Gorenstein, Shapiro & Falco 1988; Witt, Mao 
& Keeton 2000; Kochanek 2002). Given a time 
delay, H scales as 2 — (3, where (3 is the profile 
index of the potential, <fi oc r 13 . Work by Koop- 
mans & Treu (2003) shows that mass profiles may 
vary from an isothermal slope by up to 10% for 
single galaxies, producing corresponding uncer- 
tainties in H . The problem is particularly seri- 
ous for four-image systems, because the images 
are all approximately the same distance from the 
centre of the lens and thus constrain the radial pro- 
file of the lensing potential poorly. On the other 
hand, for CLASS B 1933+503, with three sources 
producing ten images, the radial mass profile is 
well constrained (Cohn et al. 2001). Unfortunately 
B 1933+503 does not show radio variability (Biggs 
et al. 2000) and is optically so faint that measuring 
a time delay is likely to be very hard. 

In some cases Einstein rings may pro- 
vide enough constraints, despite the necessity to 
model the extended source which produces them 



(Kochanek, Keeton & McLeod 2001) although 
models constrained by rings may still be degen- 
erate in H (Saha & Williams 2001). 

The mass profile degeneracy is particularly 
sharply illustrated by "non-parametric" modelling 
of lens galaxies (Williams & Saha 2000; Saha 
& Williams 2001). Such models assume only 
basic physical constraints on the galaxy mass 
profile, such as a monotonic decrease in sur- 
face density with radius. They find consistency 
with the observed image data for a wide range 
of galaxy mass models, which are themselves 
consistent with a wide range of H . Combin- 
ing two well-constrained cases of lenses with a 
measured time delay, CLASS B 1608+656 and 
PG 1 1 15+080, Williams & Saha (2000) find H = 
61 ± 18kms -1 Mpc -1 (90% confidence). 

There are a number of approaches to the res- 
olution of the mass profile degeneracy problem. 
One is to assume that galaxies have approximately 
isothermal mass distributions. (J3 ~ 1). There are 
two parts to the lensing argument in favour of 
isothermal galaxies: from the lack of odd images 
near the centre of observed lens systems Rusin & 
Ma (2001) are able to reject the hypothesis that 
significant number of lensing galaxies have pro- 
files which are much shallower {(3 > 1.2) than 
isothermal, assuming a single power-law model is 
appropriate for the mass contained interior to the 
lensed images. Similarly it can be shown that mod- 
els which are significantly steeper than isother- 
mal are unable to reproduce constraints from posi- 
tions and fluxes in existing lenses with large num- 
bers of constraints (e.g. Munoz, Kochanek & Kee- 
ton 2001; Cohn et al. 2001). The most straight- 
forward approach, that of assuming an isothermal 
lens, has been taken by many authors. In most 
cases this yields H estimates of between 55 and 
70 kms^Mpc" 1 (e.g. Biggs et al. 1999; Koop- 
mans & Fassnacht 1999; Koopmans et al. 2000; 
Fassnacht et al. 1999) but studies of some lenses 
imply much lower values (e.g. Schechter et al. 
1997; Barkana 1997; Kochanek 2003). In fact, 
Kochanek (2003) finds a serious discrepancy with 
the //ST key project value of // =71km s _1 Mpc -1 
(Mould et al. 2000; Freedman et al. 2001) unless, 
far from being isothermal, galaxy mass profiles 
follow the light distribution. 

Falco, Gorenstein & Shapiro (1985) pointed 
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out the second important problem. Any nearby 
cluster produces a contribution to the lensing po- 
tential in the form of a convergence which is 
highly degenerate with the overall scale of the 
lensing system and hence with H . Unfortunately, 
the systems with the most accurately determined 
time delays and the best-known galaxy positions 
are often those with large angular separation such 
as 0957+561, and these are the systems in which 
lensing is most likely to be assisted by a clus- 
ter. Again progress can be made by appropriate 
modelling of the cluster, and many attempts have 
been made to do this for 0957+561 (e.g. Kundic 
et al. 1997; Bernstein & Fischer 1999; Barkana 
et al. 1999) although there remain uncertainties 
in the final H estimate. As an alternative, the 
optical/infra-red images of the host galaxy may 
make an important contribution towards the break- 
ing of degeneracies (Keeton et al. 2000). 

Kochanek & Schechter (2004) summarise the 
contribution of lensing to the H debate so far 
and present options for further progress. One ap- 
proach is simply to accumulate more H determi- 
nations and rely on statistical arguments to iron 
out the peculiarities which affect each individual 
lens system; this approach is vulnerable only to a 
systematically incorrect understanding of galaxy 
mass profiles. The alternative approach is to select 
a lens system in which additional observational ef- 
fort is most capable of decreasing the systematic 
errors on the H estimate to acceptable levels. In 
this paper we argue that CLASS B0218+357 is the 
best candidate for this process. We describe new 
Hubble Space Telescope (HST) observations using 
the Advanced Camera for Surveys (ACS) which 
are aimed at removing the last major source of 
systematic uncertainty in this system. We then de- 
scribe how we use the imaging data to derive a 
value for the Hubble constant. 



2 CLASS B0218+357 AS A KEY OBJECT IN H 
DETERMINATION 

CLASS B02 18+357 was discovered during the 
early phase of the CLASS survey, the Jodrell 
Bank-VLA Astrometric Survey (JVAS; Patnaik et 
al. 1992). B0218+357 consists of two images (A 
and B) of a background flat-spectrum radio source 
separated by / /334, together with an Einstein ring 



(Patnaik et al. 1993). The optical spectrum shows 
a red continuum source superimposed on a galaxy 
spectrum. The redshift of the lensing galaxy has 
been measured optically by Browne et al. (1993) 
and Stickel & Kuehr (1993), and by Carilli, Ru- 
pen & Yanny (1993) at radio wavelengths giving 
the most accurate result of 0.6847. Cohen et al. 
(2003) have measured the source redshift of 0.944. 

It quickly became apparent that the lensing 
was performed by a spiral galaxy. Spiral lenses 
generally produce smaller image separations than 
elliptical lenses because of their lower mass. The 
spiral nature of B02 18+357 was deduced directly 
from early high-resolution optical images from 
the Nordic Optical Telescope (Grundahl & Hjorth 
1995), and was consistent with evidence from 
molecular line studies which revealed absorption 
of the radio emission from the background radio 
source by species in the lensing galaxy includ- 
ing CO, HCO+, HCN (Wiklind & Combes 1995), 
formaldehyde (Menten & Reid 1996) and water 
(Combes & Wiklind 1997). Moreover, in the opti- 
cal, the A image, which is further from the galaxy, 
is fainter than the B image (Grundahl & Hjorth 
1995) , despite being a factor ~3 brighter in the ra- 
dio. This suggests that the line of sight to the A im- 
age intercepts a great deal of dust, possibly asso- 
ciated with a giant molecular cloud in the galaxy. 
The lensing galaxy appears close to face-on, a con- 
clusion deduced from its symmetrical appearance 
in optical images. This conclusion is consistent 
with the small velocity line-width of the absorp- 
tion lines (e.g. Wiklind & Combes 1995). 

Further radio imaging resolved both the A and 
B images into core-jet structures (Patnaik, Porcas 
& Browne 1995; Kemball et al. 2001; Biggs et al. 
2003), as well as revealing more details of the Ein- 
stein ring (Biggs et al. 2001). The combined con- 
straints from the core-jet structure and the ring to- 
gether strongly constrain mass models. A and B 
lie at different distances from the galaxy, and with 
the Einstein ring constraints permits determination 
of both the angular structure of the lensing mass 
(Wucknitz et al. 2004) and (most importantly) the 
mass-radius relation for the lens. 

Biggs et al. (1999) have determined a time 
delay of 10.5±0.4 days (95% confidence) for 
B021 8+357 using radio monitoring observations 
made with the VLA at both 8.4 GHz and 15 GHz. 
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Consistent results were obtained from the varia- 
tions in the total intensity, the percentage polar- 
ization and the polarization position angle. Biggs 
et al. used the time delay and existing lens model 
to deduce a value for the Hubble constant of 
69t\l kms- 1 Mpc- 1 (95% confidence). It should 
be noted, however, that the error bars on the as- 
sumed position for the lensing galaxy with respect 
to the lensed images were over-optimistic and 
hence their quoted error on H is too small. Cohen 
et al. (2000) also observed B02 18+357 with the 
VLA, and measured a value for the time delay of 
lO.llJ g days. This corresponds to an H value of 
71I21 kms -1 Mpc~ 1 (95% confidence), the larger 
error bars in Cohen et al.'s measurement being due 
to their use of a more general model for the source 
variability, although they used the same model for 
the lensing effect as Biggs et al. The error bars do 
not take into account any systematic error associ- 
ated with the uncertain galaxy position. 

Lehar et al. (2000) used the then existing con- 
straints to model the CLASS B02 18+357 sys- 
tem. They found that, even for isothermal mod- 
els, the implied value of H was degenerate with 
the position of the centre of the lensing galaxy, 
with a change of about 0.7 km s _1 Mpc _1 (about 
1 per cent) in the value of H for every 1 mas shift 
in the central galaxy position. Their uncertainty on 
the position derived from HST infra-red observa- 
tions using the NICMOS camera is approximately 
±30 mas. 

Recently, using a modified version of the 
LENSCLEAN algorithm (Kochanek & Narayan 
1992; Ellithorpe et al. 1996; Wucknitz 2004), 
Wucknitz et al. (2004) have been able to constrain 
the lens position from radio data of the Einstein 
ring. With the time delay from Biggs et al. (1999) 
of (10.5 ± 0.4) days, they obtain for isothermal 
models a value of H = 78±6 km s" 1 Mpc- 1 (95% 
confidence). They use VLBI results from other au- 
thors (Patnaik et al. 1995, Kemball et al. 2001) as 
well as their own data (Biggs et al. 2003) to con- 
strain the radial profile from the image substruc- 
ture and obtain (3 ~ 1.04, very close to isothermal 
{(3=1). Our aim in this paper is to use new optical 
observations to determine the lensing galaxy posi- 
tion directly and to compare this with the indirect 
determination of Wucknitz et al. (2004). 

We briefly summarise the reasons why, 



given the observations presented here, CLASS 
B02 18+357 offers the prospect of the most unbi- 
ased and accurate estimate of H to date. 

(i) The observational constraints are arguably 
the best available for any lens system with a mea- 
sured time-delay. 

(ii) The radio source is relatively bright (a few 
hundred mJy at GHz frequencies) and variable 
at radio frequencies, so time delay monitoring is 
relatively straightforward and gives a small error 
(Biggs et al. 1999) which can be improved with 
future observations. 

(iii) The system is at relatively low redshift. 
This means that the derived values for H will not 
depend on the matter density parameter and cos- 
mological constant by more than a few percent. 

(iv) The lens is an isolated single galaxy and 
there are no field galaxies nearby to contribute to 
the lensing potential (Lehar et al., 2000). 

Although most lenses have at least one of these 
desirable properties, CLASS B02 18+357 is the 
only one known so far which has all of them. It 
thus becomes a key object for H determination. 
It is only the lack of an accurate galaxy position 
that led in the past to it being excluded from con- 
sideration by many authors (e.g. Schechter 2001; 
Kochanek 2003). 

3 THE ACS OBSERVATIONS 

Resolving the lens galaxy and lensed images is an 
aim that benefits from high resolution combined 
with high dynamic range, and so we asked for and 
were awarded time on the Advanced Camera for 
Surveys (ACS; Clampin et al. 2000) on the Hubble 
Space Telescope. 

Although observations in the blue end of the 
spectrum would have maximised angular resolu- 
tion, the likely morphological type of the lensing 
galaxy (Sa/Sab) meant that asymmetry due to star 
formation in the spiral arms could have caused 
problems in the deconvolution process which re- 
lies on symmetry in the lensing galaxy (see Sec- 
tion 5). Hence observations at wavelengths longer 
than 4000A in the rest frame of the galaxy were 
desirable. At the redshift of the lens this dictated 
the use of /band, i.e. the F814W filter on the HST. 

The ACS has two optical/near-IR "channels", 
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X offset (pixels) 

Figure 1. 4+4 dither pattern used for the observations of 0218+357. The 
pattern provides steps to the level of i -pixel. 

the Wide Field Channel (WFC) and the High 
Resolution Channel (HRC). The HRC's pixel re- 
sponse exhibits a diffuse halo at longer wave- 
lengths due to scattering within the CCD. As a re- 
sult, roughly 10% of the flux from a point source 
will be scattered into this halo at 8000A, possi- 
bly making it more difficult to resolve the lens- 
ing galaxy from the lensed images. We selected 
the WFC for use in our observing programme 
since it does not suffer from this effect. Unlike 
the HRC, the WFC moderately under-samples the 
HST point-spread function (PSF) at 8000A. To 
counter this effect, we selected two distinct four- 
point dither patterns, alternating between them 
over the course of the observations. The dither pat- 
terns used are shown in Figure [T] (see also Mutch- 
ler & Cox 2001 for more information on HST 
dither patterns). 

The WFC has a field of view of 202" x 202", 
and a plate scale of approximately 50 mas pixel -1 . 
Since we used a gain of unity, saturation in the 
images occurs near the 16-bit analogue to digital 
conversion limit of 65,000 e~ pixel -1 rather than 
at the WFC's full well point of 85,000 e pixel -1 . 
We determined the exposure time required on 
B02 18+357 through simulation. 

The full programme of B02 18+357 observa- 
tions was carried out over the period from Au- 
gust 2002 to March 2003. Details of the observ- 
ing dates and exposure times are shown in Table[T] 
The total available telescope time was split into 7 
visits on B02 18+357, 6 of which provided 2 hours 
integration time on the science target. The remain- 



ing visit (16) was designed to permit the pro- 
gramme to be salvaged in the unlikely case that the 
observing pattern chosen for visits 10-15 proved 
to be both inappropriate and uncorrectable. This 
visit suffered from increased observing overheads 
relative to the other visits and provided an integra- 
tion time of 1 hour, 22 minutes on B02 18+357 as 
a result. In order to deconvolve the images we re- 
quired an ACS/WFC PSF, so two short visits (1 
and 2 in Table[TJ were dedicated to observing two 
Landolt standard stars (Landolt 1992). Follow- 
ing McLure et al. (1999), observations were taken 
with several different exposure times to allow the 
construction of a composite PSF that would have 
good signal-to-noise in both core and wings whilst 
avoiding saturation of the core. Standards were se- 
lected to be faint enough not to saturate the WFC 
chip on short integration times, and to have the 
same V — I colour, to within 0.2 magnitudes, as 
the lensed images in the B02 18+357 system. The 
exposure times on the standard stars ranged from 
0.5 seconds to 100 seconds each, the longest ex- 
posures each being split into two 50 second expo- 
sures to simplify cosmic ray rejection. 

4 REDUCTION OF THE ACS DATA 

The uncalibrated data produced by automatic pro- 
cessing of raw HST telemetry files by the OPUS 
pipeline at the Space Telescope Science Institute 
(STScI) were retrieved along with flat fields, su- 
perdarks and other calibration files. 

The data were processed through the ACS cali- 
bration pipeline, CALACS (Pavlovsky et al. 2002), 
which runs under NOAO's IRAF software. The 
CALACS pipeline de-biased, dark- subtracted and 
flat-fielded the data, producing a series of cali- 
brated exposures. The pipeline also combined the 
CR-SPLIT exposures in visits 1, 2 and 16 to elim- 
inate cosmic rays. The calibrated exposures were 
in general of acceptable quality for use in the next 
stage of reduction, except for visit 15 in which 
there was some contamination of the images by 
stray light, probably from a WFPC2 calibration 
lamp (R. Gilliland and M. Sirianni, private com- 
munication). It is possible that this defect can be 
corrected in the future using the techniques which 
were used by Williams et al. (1996) to remove 
stray light from some HDF exposures, but we 
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Visit no. 


Target 


Observation date 


Exposure time 


Dither pattern 


File name r 


10 


CLASS B021 8+357 


2003 February 28 


20x360 sec 


4+4 


j8d410 


11 


CLASS B021 8+357 


2003 March 01 


20x360 sec 


4+4 


j8d411 


12 


CLASS B0218+357 


2003 January 17-18 


20x360 sec 


4+4 


j8d412 


13 


CLASS B021 8+357 


2003 March 06 


20x360 sec 


4+4 


j8d413 


15 


CLASS B021 8+357 


2003 March 11 


20x360 sec 


4+4 


j8d415 


14 


CLASS B021 8+357 


2002 October 26-27 


20x360 sec 


4+4 


j8d414 


16 


CLASS B0218+357 


2002 September 1 1 


12x360 sec 


4+4 


j8d416 








8x75 sec 


4C 




1 


92 245 


2002 October 17-18 


8x0.5 sec 


4+4 


j8d401 








8 x 8 sec 


4+4 










1x360 sec 






2 


PG0231+051B 


2002 August 25 


8x0.5 sec 


4+4 


j8d402 








8 x 8 sec 


4+4 










8 x 50 sec 


4C 





Table 1. Log of HST observations. All observations were taken with the ACS using the F814W filter, corresponding to / band. A 4+4 dither pattern refers to 
an eight-point dither consisting of two nested parallelograms, whereas 4C refers to a four-point dither parallelogram with the exposure at each point split into 
two for explicit cosmic-ray rejection. 



have not attempted to deconvolve the contami- 
nated visit. 

The calibrated exposures were fed to the next 
stage of reduction, based around the dither pack- 
age (Fruchter & Hook 2002), and the STSDAS 
packages pydrizzle (Hack 2002) and multidrizzle 
(Koekemoer et al. 2002). These tools clean cos- 
mic rays, remove the ACS geometric distortion 
and "drizzle" the data on to a common output im- 
age (Fruchter & Hook 2002). The drizzle method 
projects the input images on to a finer grid of out- 
put pixels. Flux from each input pixel is distributed 
to output pixels according to the degree of overlap 
between the input pixel and each output pixel. To 
successfully combine dithered images into a single 
output image, knowledge of the pointing offsets 
between exposures is needed. The expected offsets 
are determined by the dither pattern used, but the 
true offsets might vary from those expected due 
to thermal effects (Mack et al. 2003) within sin- 
gle visits. To determine the true pointing offsets 
between dithered exposures, we cross -correlated 
the images. Since the WFC at I-band moderately 
undersamples the HST's PSF, and since most of 
the features detected in the images were extended 
rather than stellar, we opted for pixel-by-pixel 
cross -correlation rather than comparisons of posi- 
tions of stars between different pointings, to max- 
imise our use of the available information. 

Images in each visit were drizzled on to a 
common distortion-corrected frame and then pairs 
of these images were cross-correlated. The two- 
dimensional cross-correlations have a Gaussian 



shape near their centres. The shift between pairs 
of images is measured by fitting a Gaussian func- 
tion to the peak in the cross-correlation, and the 
estimated random error in the shift is derived from 
the position error given by the Gaussian fit. For 
our data, the random errors in the measured shifts 
ranged from 0.8 to 2.5 mas. The RMS scatter be- 
tween corresponding pointings within a visit was 
typically less than 10 mas, or 20% of a single 
WFC pixel. We fed these shifts to the multidriz- 
zle script, which carried out the drizzling of visits 
to common, undistorted output frames. As part of 
the drizzling process, we opted to decrease the out- 
put pixel size from 50 mas square (the natural size 
of the undistorted output pixels) to 25 mas square. 
To avoid blurring and "holes" in the output im- 
age, the input pixels were shrunk to 70% of their 
nominal size before being drizzled on to the output 
frames (Fruchter & Hook 2002). We used a gaus- 
sian drizzling kernel to slightly improve resolu- 
tion and reduce blurring. At the end of this process 
each visit, except for visit 15, provided us with a 
output image with improved sampling compared 
to the individual input exposures. The deconvolu- 
tion of these images is described in Section |3 

Since deconvolution depends greatly on the 
accuracy of the PSF model, we have produced a 
number of different PSFs. Unfortunately the Lan- 
dolt standard star (Landolt 92 245) observed in 
visit 1 was resolved into a 0.5" double by the 
ACS/WFC, so we have concentrated on extracting 
a PSF from visit 2 (Landolt PG0231+051B). The 
calibrated exposures were combined using mul- 
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tidrizzle. Saturated pixels were masked. The re- 
sulting PSF suffered from serious artifacts consist- 
ing of extended wings approximately 80 mas up 
and down the chip from the central peak, possibly 
due to imperfect removal and combination of pix- 
els which were affected by bleeding of saturated 
columns. We believe these extended wings to be 
artifacts since they rotate with the telescope rather 
than the sky, and are not present in stars in other 
visits. 

In lieu of ideal standard star PSFs, we have 
created per-visit PSFs by averaging field stars to- 
gether. These field star PSFs do not suffer from the 
artifacts present in the visit 2 PSF. Manual exam- 
ination, pixel-by-pixel, of the difference between 
the PSFs and the central regions close to image 
B, indicate that the RMS error in the PSFs is be- 
tween 5 and 15%. Such an error increases linearly 
with the counts, rather than with their square root; 
lacking a perfectly-fitting PSF, we have allowed 
for this error when performing the data analysis 
described in Section 



5 ANALYSIS OF THE POSITION DATA AND RESULTS 
5.1 General remarks 

Figure El shows the ACS image of the CLASS 
B02 18+357 system from the combined dataset, 
consisting of all science visits on B02 18+357 ex- 
cluding visit 15. To produce this combined im- 
age, separate visits were related through restricted 
linear transformations (rotation and translation) 
based on the positions of unsaturated stars com- 
mon to all images, and then co-added. In locating 
the lensing galaxy, however, we did not use this 
combined image, preferring to work with the sepa- 
rate visits. The two compact images (A and B) can 
be distinguished as can the lensing galaxy which 
lies close to B. The spiral arms of the galaxy are 
clearly seen confirming the earlier deductions that 
the lensing galaxy is a spiral (Wiklind & Combes 
1995; Carilli et al. 1993). The spiral arms appear 
to be smooth and regular and there is no sign of 
significant clumping associated with large-scale 
star formation. The galaxy appears almost exactly 
face-on. We deduce this by assuming a galaxy po- 
sition close to B and comparing counts between 
pixels at 90 degree angles from each other about 




Figure 2. Combined ACS image of B0218+357. The lensed images are both 
visible; the brighter image, B, is close to the centre of the lensing galaxy. 
The spiral arms of the galaxy are clearly visible. The plot above the image 
shows a one dimensional slice passing through images A and B. The best-fit 
positions of A and B on this slice are marked, along with A' , the position 
of A expected from the radio image separation (334 mas). The separation 
between A and B in the optical image is 3 17+4 mas (2 a). 

the assumed centre. Examination of the residuals 
reveals no sign of ellipticity. 

The core of the lensing galaxy is strongly 
blended with B (Figure and is relatively weak. 
Extrapolation of an exponential disk fit to the outer 
isophotes of a slice through the central region 
shows that the peak surface brightness of image 
B exceeds that of the galaxy by a factor of about 
30-50. Thus the determination of an accurate po- 
sition for the galaxy is a challenging task. Before 
discussing the process in more detail below, we 
outline the various steps we go through to obtain a 
galaxy position. They are: 

(i) For each visit we measure the positions of 
the A and B images. 

(ii) We subtract PSFs from these positions. 

(iii) Using PSF-subtracted data we look for the 
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galaxy position about which the residuals left af- 
ter PSF subtraction appear most symmetric. We do 
not subtract a galaxy model from the images. This 
approach finds the centre of the most symmetric 
galaxy consistent with the data. 

(iv) We compute the mean and variance of the 
galaxy positions found from the individual visits. 



5.2 Analysis procedure 

Our reduced ACS/WFC images are over 
8450x8500 pixels in size, and cover 202"x202"on 
the sky. We cut out a region of 128x128 pixels 
(3.2"x3.2") centred on the lens and analyse this 
to make fitting computationally practical and to 
isolate the lens from other objects on the sky. 

Since the images, particularly image B, have 
much higher surface brightnesses than the galaxy, 
their positions can be located relatively accurately 
by subtracting parametric or empirical PSFs from 
the data and minimising the residuals. Fits were 
carried out on circular regions centred on the 
brightest pixel of each compact image. The re- 
gions chosen were 1 1 pixels in diameter. To avoid 
any bias arising from the choice of PSF, we used 
both parametric models (Airy and Gaussian func- 
tions) and the field star PSFs in the fits. The field 
star PSFs were consistently better fits to both A 
and B than the parametric models, Gaussians be- 
ing insufficiently peaked and Airy functions hav- 
ing diffraction rings that were too prominent. A 
linear sloping background was modelled along 
with the PSF in order to take account of the flux 
due to the galaxy. Typically the various methods 
agreed on positions to within a tenth of a pixel 
(2.5 mas). 

The separation of A and B determined by op- 
tical PSF fitting is consistently less than the radio 
separation of 334 mas. We find that the mean im- 
age separation in the optical is 317±2 mas (1 a) 
when the field star PSFs are used, and 315±4 mas 
when Gaussian PSFs are used. The correspond- 
ing result for Airy function PSFs is 31 1±10 mas. 
These values are mean separations taken over 
the six processed visits on B02 18+357. We have 
checked the plate scale of drizzled images against 
stars listed in the US Naval Observatory's B1.0 
catalogue, and find that the nominal drizzled plate 



scale of 25 mas pixel -1 is correct to better than 1 % 
for all visit images. 

An anomalous optical image separation has 
been suggested before for B02 18+357, starting 
with ground-based optical imaging by Grundahl & 
Hjorth (1995) and again by Hjorth (1997). Jack- 
son, Xanthopoulos & Browne (2000) used NIC- 
MOS imaging to find an image separation of 
318±5 mas, in agreement with our result from 
field star PSFs. 

We hypothesise that that this low separation 
may be a result of the high, and possibly spatially 
variable, extinction in the region of A. We suggest 
that some of the image A optical emission arises 
from the host galaxy rather than from the AGN 
which dominates the B image emission. Thus the 
centroid of A may not be coincident with the AGN 
image. Image A may be obscured completely and 
the emission seen could be due to a large region of 
star formation associated with an obscuring giant 
molecular cloud. In view of this possibility and the 
fact that B is consistently much brighter than A in 
the optical, we measure galaxy positions as offsets 
from our measured position for B. Thus, although 
A's position may be distorted in the optical it does 
not directly influence our measurements of H . 

Having determined positions and fluxes for A 
and B, we created an image containing two PSFs 
as a model for the flux from the lensed images 
alone (the "model image"). A model image was 
made separately for each visit and subtracted from 
each observed image, leaving a residual image 
which contained only the galaxy plus subtraction 
errors. Figure |3] shows a typical image with A and 
B subtracted. The model allowed us to keep track 
of how much PSF flux was removed from each 
pixel in producing the residual images. 

We opted to use the criterion of maximum 
symmetry in the residuals as a goodness-of-fit pa- 
rameter rather than attempting to fit a parametric 
model, such as an exponential disk profile, to the 
light distribution of the galaxy. The symmetry fit 
statistic is expressed as 



X 



sir — sir 



air) 2 + a(r 



A2 ' 



(1) 



in which s(r) is the count-rate at image pixel 
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Figure 3. A visit image with A and B subtracted after fitting fluxes and posi- 
tions. The residuals near the centre of each subtracted image show maxima 
of approximately 20% of the unsubtracted light. 



position r = (x, y) in counts per second, er(r) is 
the estimated noise at r, P is the set of pixels in- 
cluded in the fit (see below) and r' is the reflection 
of r around the galaxy position g, given by simple 
geometric considerations as 



r' = 2g — r. 



(2) 



The random error for each pixel in the image, 
<r(r) is estimated from the CCD equation (Merline 
& Howell 1995) summed with a contribution from 
the assumed random error in the PSF. The estimate 
of the noise is given by 



-(if = 



q 4 q- 2 Ts{v) + S + R 2 
^ N 



+ » 2 s p (r) 2 , (3) 



where r is the integration time at a single 
dithered pointing, /i is the assumed fractional er- 
ror in the PSF, s p (r) is the count-rate from A and 
B alone (stored in the model image), q is the driz- 
zling scale factor (the ratio between output and in- 
put pixel size, 0.5 for these images) N is the num- 
ber of dithered pointings combined in the drizzle 
process, S is the sky noise and R is the ACS/WFC 
read noise. Both sky noise and read noise are ex- 
pressed in units of electrons, and the integration 
time is in seconds. The resulting noise figure has 
units of counts per second, as does the drizzled im- 



age. We have checked this noise estimate against 
the background noise in our images and against 
simulations of the drizzling process. 

The set of pixels (P) included in the calculation 
of this x 2 figure can bias the fit if it is ill-chosen. 
When r' falls outside the boundaries of the image, 
the pixel r is considered to contribute nothing to 
the x 2 statistic and the pixel is not included in the 
set P. Such pixels therefore do not contribute to the 
number of constraints available and as a result do 
not increase the number of degrees of freedom in 
the fit. Alternative treatments can introduce bias; 
for instance, if these pixels are assigned large x 2 
values the fitting program is biased towards plac- 
ing the galaxy in the geometrical centre of the im- 
age. If the same pixels are considered to contribute 
zero towards the x 2 statistic but are still counted 
as part of the set P, they increase the number of 
degrees of freedom in the fit and bias the fit to 
positions away from the image's geometrical cen- 
tre. To avoid these possibilities we do not count 
degrees of freedom from pixels whose reflection 
about the galaxy centre ends up outside the image 
boundaries. 

The symmetry criterion is non-parametric and 
has the advantage of minimizing the assumptions 
that are imposed on the data; the use of a particu- 
lar distribution as a function of radius in any case 
contains an implicit assumption of symmetry. Us- 
ing the symmetry criterion on its own is in princi- 
ple robust whether or not the galaxy has a central 
bulge, and should also be unaffected if the galaxy 
contains a bar. The symmetry criterion will also 
hold for galaxies having moderate inclinations to 
the line of sight. For a circularly symmetric galaxy, 
the main effect of a small deviation away from a 
face-on orientation will be to render the observed 
image slightly elliptical. The basic symmetry cri- 
terion is that points and their reflections about 
the true centre of the galaxy's image should have 
the same flux (to within the measurement errors). 
Therefore whether the galaxy's image has circu- 
lar or slightly elliptical isophotes is unimportant 
because in both cases the same isophote passes 
through both point and reflection. This argument 
breaks down for spirals with significant inclina- 
tions as absorption is likely to become important 
and destroy any symmetry present in the image. 

The centre of maximum symmetry could be 
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displaced by spiral arms if they are not themselves 
symmetric about the galaxy centre. To analyse the 
possible effect on H of spiral arms we have made 
two sets of fits. In the first set we used all the data 
and did not apply any masking. In the second set 
we masked off the most obvious spiral arms by 
using an annular mask centred on image B with 
an inner radius of 0.375"and an outer radius of 
0.875". Regions within the inner radius or beyond 
the outer radius of the mask were left free to con- 
tribute to the symmetry fit. We implemented the 
symmetry fits so that masked pixels did not con- 
tribute to either the number of degrees of freedom 
or to the x 2 value. 

The PSF error (jj, in equation |3j) can cause a 
systematic change in galaxy position when varied 
between 0.05 and 0.15 (5 to 15%). An increase 
of fi from 0.05 to 0.15 can increase H by up to 
10 km s _1 Mpc -1 . For values above 0.15 the sys- 
tematic change in galaxy position is small com- 
pared with the random error. We estimate the PSF 
error separately for each visit, by taking the range 
of the highest and lowest residuals and dividing 
that range by the peak count-rate of image B. We 
find values between 0.07 (visit 11) and 0.19 (visit 
10) for //. For the other visits, the estimated value 
for jx is found to be 0.12. In the remainder of the 
paper we do not allow [i to vary freely but fix it to 
these estimated values. 

5.3 Extraction of the galaxy position 

In applying the symmetry criterion to 0218 we 
calculated the symmetry x 2 statistic (i.e. that of 
equation [l]) for a grid of galaxy positions extend- 
ing 20 mas east to 100 mas west of B, and from 
80 mas south to 50 mas north of B. The spacing 
between adjacent grid points is 5 mas. We present 
these grids in Figure |4] for visits 10, 11, 12, 13, 14 
and 16. The grids are shown both with and without 
masking of spiral arms. Table 2 lists the resulting 
galaxy positions. 

In Figure |5] we show the per-pixel x 2 contri- 
butions between the two cases of no masking and 
masked spiral arms for the optimum galaxy po- 
sition, together with the contributions when no 
masking is used and zero PSF error is applied. 
With zero PSF error it is clear that the residu- 
als from the subtraction of A and B dominate the 



Visit Centre Centre 





(No masking) 


(Spiral arms 


masked) 




Aa 


Ad 


Aa 


A<5 


10 


+50 


+6 


+70 


+ 12 


11 


+60 


-4 


+69 


-18 


12 


+59 


+9 


+84 


+8 


13 


+54 


-2 


+72 


-5 


14 


+59 


+0 


+76 


-16 


16 


+61 


-6 


+79 


-14 


Mean 


+57+4 


+ 1+6 


+75+6 


-6+13 



Table 2. Derived optical centre of the galaxy, expressed as offsets in mas 
from the measured optical position of B. RA offsets are given with west as 
positive. 
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Figure 5. These images show the contribution of each pixel to the symme- 
try x 2 f° r visit 11, using the best-fit optical galaxy positions. The top-left 
image shows the per-pixel \ 2 contributions when no masking is used and 
no PSF error is used. The image on the top-right shows the contributions 
when no masking is applied and a 7% PSF error is used. The bottom-left 
image shows the contributions when spiral arms are masked, with no PSF 
eiTor. The bottom-right image shows the \ 2 contributions when the promi- 
nent spiral arms are masked and a 7% PSF error is used. All images are 128 
pixels (3.2") in both width and height. 

X 2 measure. The effect of a non-zero PSF error 
is to suppress the A/B subtraction residuals and 
cause the spiral arms to dominate the x 2 mea- 
sure, unless they are masked. Because it is unclear 
which position (masked or unmasked) best repre- 
sents the mass centre of the lensing galaxy, we re- 
port both masked and unmasked galaxy positions 
(and hence estimates for H ) on equal terms. 

Deriving errors on position from the individ- 
ual visits is difficult because the symmetry x 2 in- 
creases very rapidly away from the minimum. An 
error measure derived from the shape of the min- 
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Figure 4. \ 2 grids for the galaxy position. The visit number is shown to the right of each pair of grids. The right-hand plot for each visit shows the effect of 
masking out the spiral arms, whilst the left-hand plot shows the \ 2 grid when no masking is applied. The position of B is marked, as is a line pointing towards 
the (radio) A component. The ellipse represents the position of the galaxy centre found by Wucknitz et al. (2004) using LENSCLEAN modelling of the Einstein 
ring, and the dotted lines represent Hq of (90,80, 70,60, 50)km s~ 1 Mpc" 1 from left to right, assuming an isothermal model. The axes are RA/Dec offsets from 
the position of B, expressed in arc-seconds. The RA offset is given with west as positive. 



imum for a single visit implies a spuriously high 
accuracy for the galaxy position. It is likely that 
the number of degrees of freedom in the fit is over- 
estimated and that many pixels do not contribute 
any useful information to the fit statistic, since the 
drizzling process introduces correlations between 
neighbouring drizzled pixels. However, the scatter 
between positions derived from different visits is 
large. We therefore estimate errors on the galaxy 
position by taking ellipses that enclose 68% and 
95% of the measurements from all visits to define 
our 1 a and 2 a confidence levels. Figure |6] shows 
the 95% confidence ellipses on the galaxy position 
for both sets of fits, as well as the position derived 
from Lens Clean applied to VLA data by Wuck- 
nitz et al. (2004). 

The positions obtained from the symmetry fits 
were combined with the extra constraints avail- 
able from the VLBI substructure described in Pat- 
naik et al. (1995) and Kemball et al. (2001), which 



were used to constrain mass models by Wucknitz 
(2004). The optical galaxy position was combined 
with the models of Wucknitz (2004) by adding x 2 
values for the galaxy position to the x 2 values from 
the lens models. However, the x 2 values taken 
from the symmetry fitting grid have too many de- 
grees of freedom, and so we assume that the opti- 
cal position minimum is parabolic and form a new 
X 2 statistic based on our 68% and 95% confidence 
ellipses. We define the new x 2 statistic to have a 
value of 2.31 on our 68% confidence ellipse, and 
a value of 5.99 on the 95% ellipse, and sum this 
statistic with that from the lens modelling. We em- 
phasize that the scatter between visits dominates 
the random error budget for our measurement of 
Ho- 

Combining the VLBI and optical constraints 
shows that the best-fit galaxy position and the opti- 
cal galaxy position are not coincident, as shown in 
Table|3] The galaxy position shifts by up to 13 mas 
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between the optical fit and the optical+VLBI fit. 
The shapes of the confidence regions are also al- 
tered. However, the value of H is not very sensi- 
tive to this mainly northerly shift. 

6 EXTRACTION OF H 

The general relation between the time delay At it j 
between the i th and j th images, the Hubble con- 
stant H and the lens model, parametrised by the 
potential ip, is given by 

a+ 1 + z d d d d s 

where z d is the redshift of the lens, d d and 
d s are the angular size distances to the lens and 
source, respectively, d ds is the angular size dis- 
tance to the source measured from the lens, and 
4>i is the scaled time delay at the position of the i th 
image (0*), 

^^iVVWr-VW ■ (5) 

The angular size distances are normalized in 
these equations, since they do not include factors 
of H . For general isothermal models without ex- 
ternal shear the relation becomes particularly sim- 
ple and can be written as a function of the image 
positions alone, without explicitly using any lens 
model parameters (Witt et al. 2000): 

0i = ^|0i-0 o | 2 (6) 

Here 6 is the position of the centre of the lens. Ex- 
ternal shear 7 changes by a factor between 1 ±7 
depending on the relative direction, typically re- 
sulting in similar factors for the value deduced for 
Ho. A general analysis for power-law models with 
external shear can be found in Wucknitz (2002). 

Using the recipe described in previous sections 
our lens position translates to a Hubble constant of 
H = 79 ±7 km s _1 Mpc _1 in the shearless isother- 
mal case 2 without masking, and to 66 ± 9 with 
masking. 

2 A concordance cosmological model with f! = 0.3 and A = 0.7 and a 
homogeneous matter distribution is used for the calculation of all distances 
in this paper 



Estimates of external shear and convergence 
from nearby field galaxies and large scale struc- 
ture are of the order 2 per cent (Lehar et al. 2000) 
and would affect the result only to the same rel- 
ative amount, sufficiently below our current error 
estimate to allow us to neglect these effects. 

The value of the Hubble constant we derive 
depends on the slope of the mass distribution of 
the lensing galaxy. In Figure |7] we show the per- 
mitted values of the Hubble constant for different 
models - isothermal and with a variable (3 in an el- 
liptical potential model - plotted against measured 
galaxy position. The elliptical power-law potential 
is parametrised as follows: 

m = ^-r!{6) , (7) 

e (1 + e) 2 + (1-e) 2 ' W 
= (6 x ,6 y ) , (9) 

where is the Einstein radius of the model, (3 is 
the power-law index of the potential's profile and 
e is the ellipticity of the potential. For details of 
our modelling procedure the reader is referred to 
Wucknitz et al. (2004). It is evident that the pre- 
ferred value of the Hubble constant is somewhat 
reduced compared to what is obtained by forcing 
the mass distribution to be isothermal. We also 
show contours of the radial power law (5 plotted 
against galaxy position. The optical lens position 
gives =1.1318$ (2 o). 

As discussed before, B02 18+357 has the ad- 
vantage of clear substructure in the two images 
which can be mapped with VLBI. The VLBI data 
can be used independently to derive the slope 
of the mass profile of the lens (Wucknitz et al., 
2004). Biggs et al. (2003) and Wucknitz et al. find 
a value of /3=1.04±0.02. Combining the VLBI 
constraints with the optical lens position gives 
a value of /3=1.05±0.03 (95% confidence). We 
therefore adopt this value for the mass profile's 
logarithmic slope and obtain a Hubble constant of 
70±5 km s _1 Mpc~ 1 (95% confidence) for the case 
with no masking and 61 ±7 kms -1 Mpc _1 (95% 
confidence) when the spiral arms are masked. The 
ellipticity of the potential is small (about 0.04) in 
each case, but the ellipticity of the mass distribu- 
tion will be about three times this, 0.12. The lens 
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Figure 6. The optical galaxy position compared to that determined by Wucknitz et al. (2004) using LENSCLEAN. The error ellipses are 95% confidence 
regions. The left-hand plot is for the case with no masking. The right-hand plot represents the case with masking of spiral arms. The position of B is marked, 
as is a line pointing towards the A component. The dotted lines are contours of Ho in the strictly isothermal case, and correspond (from right to left) to Ho = 
(90,80,70,60,50) kms -1 Mpc -1 . The positions determined from optical data are shown as bold ellipses. 



Data used 



Optical 

Optical 

VLBI+Optical 

VLBI+Optical 

VLBI+Optical 

VLBI+Optical 



Masking Mass profile Position (mas) Ho Ho 

Aa AS (isothermal) (variable (3) 



None 
Spiral 
None 
Spiral 
None 
Spiral 



Isothermal 
Isothermal 
Variable 
Variable 



+57 
+75 
+60 
+74 
+60 
+74 



+0 

-5 
-13 
-19 
-12 
-18 



79±7 
66+9 
74+5 
64+7 



68+6 
56 



+ 12 

15 



70+5 
61+7 



, ,0+0.07 
1 - 1J -0.09 

1.16+0.19 



1.05+0.03 
1.05+0.04 



Ellipticity 
(of potential) 

0.08+0.03 
0.05+0.04 
0.05+0.02 
0.03+0.02 
0.04+0.02 
0.04+0.02 



Table 3. Lens galaxy positions and the corresponding values of Ho and the mass profile slope (3. The optical positions are derived from the ACS images only. 
The "VLBI+Optical" positions incorporate constraints from the LENSCLEAN-based lens modelling of Wucknitz et al. (2004). The "Mass profile" column 
indicates what mass profile was assumed when combining the VLBI and optical constraints. Ho values are given in kms -1 Mpc -1 . Position offsets are 
referenced to image B, and RA offsets are given taking west as positive. All errors are quoted at 95% confidence. 



galaxy could be more inclined than it appears, or 
there could be a bar or similar feature present. 



7 CONCLUSIONS 

We have analysed the deepest optical image yet 
taken of B02 18+357 to measure the position of the 
lens galaxy. We find that simple subtraction of a 
parametric galaxy model and two point sources is 
insufficient to constrain the galaxy position, and 
we confirm earlier suggestions that the image sep- 
aration in the optical is lower than that in the radio, 
most probably due to significant extinction around 
image A. Taking advantage of the symmetric ap- 
pearance of the lens, we have defined the centre 
as that point about which the residuals (after sub- 
traction of A and B) are most symmetric. To ac- 
count for artifacts in our empirical PSF model we 
have introduced an extra noise term. We have also 
masked off the most prominent spiral arms to test 



the effect on H . We find that the lens galaxy po- 
sition is 57 ±4 mas west and 0±6 mas south of 
image B when no masking is applied. Combined 
with the results of Wucknitz et al. (2004) this leads 
to a value for H of 70±5 km s" 1 Mpc- 1 (95% 
confidence). When the most obvious spiral arms 
are masked out, we find an optical galaxy posi- 
tion of 75 ±6 mas west and — 5±13 mas south 
from image B. This results in a value for H of 
61 ±7 km s _1 Mpc~ 1 (95% confidence) when com- 
bined with VLBI constraints. 



Further work on this lens will involve in- 
creased use of LensClean to further limit the 
power law exponent (3 using VLBI constraints. 
Observations have also been made using the VLA 
with the Pie Town VLBA antenna, which together 
with VLBI will further improve the lens model for 
this system. 
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Figure 7. The 95% confidence regions for the cases of no masking (top row) and masked spiral arms (bottom row). The plots in the left-hand column show 
the galaxy position confidence regions and contours of Ho calculated for isothermal models. The plots of the central column show the confidence regions 
superimposed over contours of Ho calculated whilst allowing the logarithmic mass slope fi to vary. This reduces Ho slightly compared to the fixed f3 = 1 
(isothermal) case. The plots in the right-hand column show the same confidence regions superimposed over contours of 0. The confidence regions take both 
our optical position and the VLBI constraints used by Wucknitz et al. (2004) into account. The 95% confidence error ellipse of Wucknitz et al. is also shown. 
Contours of Ho are again at (90,80,70,60,50)kms _1 Mpc -1 reading from left to right across an image, and contours of are (0.8,0.9,1.0,1.1,1.2,1.3) reading 
from bottom to top. The optical+VLBI error ellipses are shown in bold relative to the LENSCLEAN ellipse of Wucknitz et al. 
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